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Abstract 

Particles with millicharge q and sub-eV mass can be produced in photon- 
photon collisions, distorting the energy spectrum of the Cosmic Microwave 
Background. We derive the conservative bound q < 10 _7 e (as well as model- 
dependent bounds two orders of magnitude stronger), incompatible with pro- 
posed interpretations of the PVLAS anomaly based on millicharged production 
or on millicharged-mediated axion-like couplings. 

1 Introduction 

New particles q with mass m q and charge gCe (thereby called 'millicharged') can naturally 
arise in field-theory models, e.g. via small kinetic mixing of the photon with a new vector. 
Distinctive aspects and constraints on millicharged particles are reviewed in [1]. 

Renewed interest in millicharged particles was prompted by an experimental result claimed 
by the PVLAS collaboration [2]: the polarization of a linearly polarized laser beam (photons 
with ~ 1 eV) rotates ('dichroism') and develops a tiny elliptical component ('birefringence') 
after multiple passages trough a vacuum Fabry- Perot cavity containing a (rotating) 5.5 Tesla 
magnetic field orthogonal to the beam direction. If confirmed, PVLAS calls for new physics, 
possibly pointing to the existence of new light particles that, interacting with photons, drain 
part of the beam energy in a polarization-dependent way. 

Thanks to the external field, photons in the beam can non-perturbatively convert into couples 
of new light millicharged particles: one can fit PVLAS for charge q ~ 3 • 10 _6 e and mass 
m q ~ 0.1 eV [3, 4]. Alternatively, the PVLAS birefringence could be due to mixing, induced 
by the magnetic field, of laser photons with a new axion-like scalar a, coupled to photons by 
an effective interaction term (a/ AM) Fa V F^ v (a pseudo-scalar axion coupled as (a/ AM) F^F^ 
produces birefringence with sign opposite to the one observed by PVLAS) [5], and the PVLAS 
dichroism could be due to conversion of some photons in axions escaping from the apparatus. 
The needed value of M ~ 2 • 10 5 GeV is excluded by a production in stars [6], unless the effective 
coupling is mediated by a loop of light enough particles, 1/M ~ aq 2 /v, that therefore need to 



have a small millicharge q ~ 10 6 e • ^yv/eV, where v is a loop function, related and comparable 
to m g [7].* 



Light millicharged particles therefore are a key ingredient of proposed new physics that can 
fit the PVLAS anomaly. Still, one needs to circumvent the too strong constraints on millicharged 
particles in [1, ?] by designing models that exploit the fact that PVLAS involves 0{ eV) energies, 
while the constraints in [1] come from physics at O(keV) energies (the temperature in the core 
of stars) or higher. An energy-dependent form factor that suppresses millicharged couplings 
at high energies can be obtained in the following way. (Although we use a basis that makes 
its presentation simpler, we are just outlining the model of [7]). The gauge group is U(l) em ® 
U(l)<g>U(l)' where the first U(l) cm is the usual electromagnetism, and new light particles (either 
scalars or fermions) are charged under the last U(l)', with charge q', assumed to be comparable 
to e [7] . The propagator matrix of the three vectors is assumed to be [7] 



i.e. canonical kinetic terms, plus a small kinetic mixing between the first two U(l), plus a vev 
that breaks the last two U(l) to their diagonal sub-group, giving a mass ~ m to their symmetric 
combination. At q 2 <C m 2 , the new light particles get a millicharge q ~ eq' under the photon. 
At q 2 3> m 2 one can neglect m 2 and this millicharge disappears: choosing m < eV allows to 
avoid the bounds from higher-energy physics discussed in [l].t 

Searches for this kind of millicharged particles have been performed using reactor neutrino 
experiments, obtaining the bound q < 10~ 5 e [10]. New reactor experiments have been discussed 
as a way to test interpretations of the PVLAS anomaly that involve millicharged particles [10]. 

We show, trying to be as model-independent as possible, that existing models are not com- 
patible with known CMB physics after recombination, which probes the same sub-eV energies 
relevant for PVLAS, and we obtain new bounds on millicharged particles. 

* Other interpretations not involving millicharged particles have been put forward. Ref. [8] suggested that 
the axion-like scalar might be replaced with a massive vector B M with a Chern-Simons-like coupling to photons 
A^. However we do not see how the <-> B M oscillations proposed in [8] can produce the dichroism claimed 
by PVLAS. It could be produced by emission of the longitudinal component of the extra vector, with couplings 
enhanced by the polarization vector e M ~ q^/m, but (just like in the axion case) this possibility is not compatible 
with bounds from star cooling. Extra observables related to axion-like physics are discussed in [9]. 

^The Higgs b! that breaks U(1)®U(1)' to its diagonal subgroup providing the vector mass m gets a millicharge 
q ~ eq' with no form factor that suppresses it at high energy: therefore the model is excluded if h' is light enough 
to be produced in stars or during BBN. Since my ~ Xm/g' , a heavy enough h' needs a quartic scalar coupling 
X\h'\ 4 with A > 1. This may appear crazy, and indeed such a non-perturbative Higgs coupling would be excluded 
if h' would break a non-abelian gauge group, because non-abelian vectors feel it (for example the SM Higgs must 
be lighter than <inMz/g ~ TeV). However the model only involves abelian vectors, and in the abelian case a big 
A remains confined to the Higgs sector, so that decoupling the physical Higgs is not impossible. Indeed this is 
well known [12]: an abelian theory remains sensible if gauge invariance is broken by adding a vector mass term 
(i.e. no Higgs or infinitely heavy Higgs). In conclusion, the models of [7, 3] an unusual but acceptable ingredient. 
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2 Model-independent cosmological constraints 

Millicharged particles affect CMB cosmology in many ways. We list the relevant processes, 
discuss their rate, their effect, and their model-independence. We restrict our attention to light 
millicharged particles, m q < eV, and on cosmology after decoupling, T< eV: according to the 
SM, photons behave as free particles with negligible thermal mass: their dispersion relation is 
q 2 = 0. We conservatively assume that the universe initially contains only ordinary SM particles. 

a) 77 — ► qq is the main process for production of millicharged particles. Being induced by a 
small adimensional gauge coupling, this process is maximally effective at low temperatures, 
~ max(m g ,To) (where To = 2.7K = 0.23 meV is the present temperature), so that its 
rate is dominated by model-independent gauge interactions, controlled only by m q and q. 

We can neglect the process — > eqq, that contains Bethe-Heitler-like contributions a ~ 
e 2 (7 4 /(4"7r) 3 X^ (suppressed with respect to 77 — ► qq by a factor ~ e 2 /(4-7r) 2 ) and Compton- 
like contributions a ~ e 4 g 2 /(4ir) 3 m 2 (suppressed by the electron mass). Both contributions give 
rates additionally suppressed by the small number density of free electrons, n e /?i 7 <C 1. 

We can estimate the number density n q of millicharged particles plus their anti-particles 
produced by process a) as 

y = ^ ~ min (V ^(77 -» qq) ff^y ) where ^(77 -» qq) ~ ( 2 ) 

and rij(T) = 2£(3)T 3 /7r 2 is the photon number density, used to normalize n q such that their 
relative density Y = n q /n^ is not affected by the expansion, with Hubble rate H(T) = R/R. 

Fig. 1 shows isocurves of Y, more precisely computed solving the Boltzmann equations 
described in the Appendix. If Y -C 1 the process 77 — > qq leads to an energy-dependent 
depletion of the CMB, whose energy spectrum has been measured with ~ 10~ 5 accuracy by 
FIRAS [11]. As illustrated in fig. 2, by fitting FIRAS data we find 

Y = n g /n 7 £ 6 • 10~ 5 at 3a confidence level {A X 2 = 9). (3) 

for m q ~ 0.1 eV, and a slightly weaker bound at smaller m q . Fig. 1 shows the precise constraint 
in the (q, m q ) plane. We ignored effects possibly related to the formation of bound states at 
recombination, as they would appear at the frequencies of the H and 4 He recombination: a region 
of the CMB spectrum not tested by FIRAS and dominated by the galactic signal of thermal 
dust. 

Before concluding that we can safely exclude the interpretation of the PVLAS anomaly based 
on [3] (the interpretation of [7] employs values of (q, m q ) that can be not far from the allowed 
region) we need to assess if this bound is model- independent, or if one can add some other new 
physics and use it to restore a thermal CMB spectrum. 

We do not see how additional interactions, slower than the expansion rate, could provide 
such a restoration since the CMB spectrum has been measured in a sizable range of energies. 
New interactions, enough faster than the interaction rate, would re-thermalize photons, but at 
the expense of thermalizing also some new particle (this happens e.g. in the model of [4]): data 
on cosmological anisotropics disfavors the resulting significant depletion of the photon energy 
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Figure 1: Isocurves of the cosmological abundance n q /n^ of fermionic (left plot) or scalar (right 
plot) millicharged particles as function of their mass m q and charge q. The shaded region is 
excluded by the CMB energy spectrum. The dot-dashed curve is the (robust, but not fully model- 
independent) constraint estimated in eq. (5), for q' ~ e. For comparison, the PVLAS anomaly 
can be interpreted as production of millicharges with q ~ fewx 10 -6 e and maybe m q ~ 0.1 eV [3]. 

density with respect to other components, and the addition of new quasi-relativistic interacting 
particles, that behave as a non freely streaming fluid. 

This problem is avoided if the new particle, added to keep photons in thermal equilibrium, 
has a mass m > Toi such that it decouples at T < m. In the simplest scenario the millicharged 
particles themselves could do the job. 

Indeed the process a) is much faster than the expansion rate if q > 10 -5 e; however such a 
large millicharge cannot fit PVLAS (that actually excludes a too large millicharge) and leads to a 
thermalized rate also at higher temperatures, thereby significantly distorting CMB anisotropies 
imprinted around last scatterings, at T ~ 0.25 eV. 

Without performing the necessary dedicated analysis, it seems unlikely that such a big modi- 
fication of cosmology could be compatible with observations: in the tight coupling limit the pho- 
ton/baryon/millicharged fluid would have a sound speed different than the usual photon/baryon 
fluid (whose acoustic oscillations have been observed as peaks in the CMB anisotropy angular 
spectrum). 

Similar problems arise if one instead assumes that millicharged particles are thermalized at 
the same temperature as photons: extra non-freely streaming particles, photons that remain cou- 
pled after recombination, a non standard amount of relativistic energy density at recombination. 
Since interaction rates mediated by adimensional couplings become slower than the expansion 
rate at large temperatures, this alternative scenario is unnatural and we see no preferred way 
that realizes it in a clean way. One can imagine many different possibilities; for example the 
effective number of neutrinos (commonly used to parameterize the relativistic energy density 
at recombination, probed by CMB data) typically deviates from 3 by some O(few) factor and 
can even be smaller than 3, in scenarios where millicharged particles reheat photons below the 
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Figure 2: FIRAS data compared to the energy- dependent depletion of the CMB spectrum due to 
77 ~~ > QQ- We plot 1 — 7*, with r given in eq. (10), computed for fermion (scalar) millicharges 
with m q = 0.1 eV and q = 10~ 7 e (q = 1.7 10~ 7 e / ) ; chosen such that the two cases give roughly 
equal effects, excluded at about 3 standard deviations. 



decoupling temperature of neutrinos. 

Once that enough millicharged particles have been produced, photons start having significant 
interactions with them. Let us now turn to examine the potentially most remarkable of such 
scattering processes. 

b) — ► 7^ leads to a photon attenuation length. Its cross section is a = q /67rm q in the 
non-relativistic Thompson limit. This process does not lead to new significant constraints 
and can be suppressed in a model-dependent way, as discussed at point d). 

Furthermore, there can be model-dependent scatterings involving new-physics particles. The 
presence of at least one new extra new light vector 7' seems almost model independent; e.g. in 
the model outlined in the introduction the millicharged particle has a sizable gauge coupling 
q' ~ e to two extra vectors 7' with masses m! = and m' ~ m. Therefore, unless all 7' are too 
heavy for being relevant, they have sizable effects. 

c) 7Q — ► -y'q, with cross section a = q 2 q' 2 /6-Krn 2 , in the limit of non-relativistic millicharge q 
and of ultra-relativistic 7'. If q' ~ e one gets a constraint orders of magnitude stronger 
than previous bounds by demanding that 75 — ► ^y'q does not distort the CMB energy 
spectrum. (This process can easily be in thermal equilibrium, again restoring a thermal 
CMB spectrum but also generating the other problems previously discussed). Coherent 
forward scattering, i.e. 7^7' oscillations, are suppressed by a mixing angle 9 ~ q/q' 
which is small in the most plausible part of the parameter space, q' ~ e. 

d) Scatterings b) and c) involve an initial state millicharge q: its density can be suppressed 
by qq — ► ■ A sizable suppression of n q is obtained only if my <T <C m q and if 
this process is in thermal equilibrium (this is indeed the most natural situation, since 
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a (QQ — *• ~ <l' 4 / m q is n °t suppressed by any millicharged coupling), while n q is only 
reduced by an 0(1) factor at T ~ m g . 

The above considerations suggest another robust (but not fully model-independent) constraint, 
coming from the distortion of the CMB energy spectrum due to "fq — ► j'q at temperatures 
T ~ T*. The resulting total deficit of photons is estimated to be 

Y 7 ~ min (l,a(jq^ ^' ^^{T) ^) Wher6 ff ^ 79 ~* ~ £rT2 ^ 
Since this effect is energy-dependent, FIRAS data imply Y 7 < 10 -4 which translates into 

1/3 / \ 1/6 



q'^W^e^-j max ^>^J ■ ( 5 ) 

This strong bound depends weakly on the model-dependent parameter q', which is comparable 
to e if one wants to get millicharged particles from a small 7/7' kinetic mixing, rather than 
putting by hand a small gauge coupling. This bound is grossly incompatible with the values of 
q, q' used in the model of [7] to mediate at one loop the effective operator aF fll/ F^ u /AM. Indeed 
this operator generates itself a Y 7 ~ 10~^ 7_8 \ a few orders of magnitude below the sensitivity of 
FIRAS, for the values of M ~ few • 10 5 GeV and of the axion-like mass m a ~ few meV suggested 
by PVLAS. 



3 Conclusions 

Motivated by the PVLAS anomaly [2], the authors of [7, 3] proposed new models containing 
millicharged particles whose charge only appears at low energy (or more precisely at small 
momentum transfer), avoiding the bounds on millicharged particles from higher-energy probes, 
such as star cooling and BBN [1]. Millicharged particles are used to directly fit PVLAS in [3], 
and to mediate an axion-like coupling to a new light scalar in [7]. 

Cosmology at low energies (namely at T eV after recombination) provides a direct and 
sensitive probe to this sort of new physics. Within standard cosmology photons behave as 
free particles with negligible thermal mass. New physics processes like 77 — > qq distort the 
CMB energy spectrum, potentially conflicting with FIRAS that measured a black-body energy 
spectrum up to ~ 10~ 5 accuracy [11]. The resulting constraints on millicharged particles, 
summarized in fig. 1, are therefore based on simple and fully safe cosmology, and exclude the 
interpretations of the PVLAS anomaly proposed in [7, 3]. We cannot propose modifications of 
the models that avoid the conflict with FIRAS data. 

Acknowledgements We thank R. Barbieri, A. Macchi, L. Maiani, E. Gabrielli, U. Gastaldi, 
M. Roncadelli and especially S. Davidson for many useful discussions and suggestions. 

Appendix 

The evolution of Y = n q /rij is described by the Boltzmann equation 
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where z = m q /T, H is the expansion rate. The thermally averaged interaction rate for 77 <-> qq is given 

by 



7A 
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ds s 1/2 Ki 
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where the reduced cross section, defined as da/dt = ^ |^4| 2 / 
final-state indices, is 



a{s) 



-/3x(4 + x) + (ar + 4x - 8) In 



a(s) (7) 
7rs where the sum runs over all initial- and 
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if millicharged particles are a Dirac fermion and 

„4 



27TX 2 



#r(4 + x) + 4(2 - x) In i±| 



(8) 



(9) 



if millicharged particles have spin 0. Here x = s/M 2 and (3 = y/l — 4/x is the q, q velocity with respect 
to their center-of-mass frame. 

Assuming F C 1 the Bolztmann equation for the distortion in the photon energy spectrum, r(e) = 
/ 7 (e)//BE(e), is 



Hz 



dr{e) 
dz 



— ^- / dede' e'f BE (e') ■ a(s = 2ee'(l - c)T 2 
62-K z e J 



(10) 



where e = E/T is the comoving photon energy, /be(c) = l/( e<! — 1) IS the Bosc-Einstcin distribution 
function and c is the cosine of the scattering angle between the two photons. (Y is computed using 
cq. (7), where, as usual, all thermal distributions are approximated with the Maxwell distribution: this 
approximation makes almost a factor of 2 difference in Y, that was compensated in eq. (3) by properly 
rescaling the right-handed side). 
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